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ABSTRACT 

We present a measurement of the galaxy merger fraction and number density from observations in 
the Hubble Ultra Deep Field for 0.5<z<2.5. We fit the combination of broadband data and slitless 
spectroscopy of 1308 galaxies with stellar population synthesis models to select merging systems based 
on stellar mass of > 10^" Mq. When correcting for mass incompleteness, the major merger fraction 
is not simply proportional to (1 + z)™, but appears to peak at Zf^ac — 1-3 ± 0.4. From this merger 
fraction, we infer that ~ 42% of massive galaxies have undergone a major merger since z ~ 1. We 
show that the major merger number density peaks at zjens — 1-2, which marks the epoch where major 
merging of massive galaxies is most prevalent. This critical redshift is comparable to the peak of the 
cosmic star formation rate density, and occurs roughly 2.6 Gyr earlier in cosmic time than the peak 
in the number density of the X-ray selected active galactic nuclei. These observations support an 
indirect evolutionary link between merging, starburst, and active galaxies. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: interaction 



1. INTRODUCTION 

The hierarchical formation scenario suggests that 
galaxy merging is the driving mechanism in galaxy for- 
mation (eg. White & Frenk 1991; Kauffmann, White, 
& Guiderdoni 1993). Since galaxy merging has been 
linked to numerous astrophysical processes, including 
star formation and the fueling of a supermassive black 
hole (SMBH), which triggers an active galactic nucleus 
(AGN), the merger history of galaxies is critical for a 
complete understanding of galaxy evolution. 

The redshift evolution of the galaxy merger fraction 
has been well studied for z < 1, and is roughly propor- 
tional to (1 -I- z)™, where typically 1 < m < 4 (Zepf & 
Koo 1989; Burkey et al. 1994; Carlberg, Pritchet, & In- 
fante 1994; Yee & EUingson 1995; Neuschaefer et al. 1997; 
Patton et al. 1997; Le Fevre et al. 2000; Patton et al. 
2002; Lin et al. 2004; Xu, Sun, & He 2004; BeU et al. 
2006; Kartaltepe et al. 2007). Toomre (1977) predicted 
the galaxy pair fraction to be proportional to the matter 
space density at low redshift, and Carlberg (1991) refined 
this idea to show that the power-law index should be 
somewhat sensitive to the Cosmological Constant. Based 
on iV-body simulations, Berrier et al. (2006) argue that 
while the galaxy merger rate per halo is increasing with 
redshift, the number of halos massive enough to host 

* BASED ON OBSERVATIONS MADE WITH THE NASA/ESA 
HUBBLE SPACE TELESCOPE, OBTAINED FROM THE DATA 
ARCHIVE AT THE SPACE TELESCOPE SCIENCE INSTL 
TUTE, WHICH IS OPERATED BY THE ASSOCIATION OF 
UNIVERSITIES FOR RESEARCH IN ASTRONOMY, INC., UN- 
DER NASA CONTRACT NAS 5-26555. 
Electronic address: russell.ryanjr@asu.edu 



a galaxy pair is decreasing. This balance will result in 
a flattening of the high redshift pair counts consistent 
with some high redshift observations (eg. Lin et al. 2004). 
Fakhouri & Ma (2007) suggest that this power-law form 
is a particular representation of a universal merger rate 
which scales as a power-law with the progenitor mass 
ratio and depends only weakly on the halo mass. Addi- 
tionally, (Conselice 2006) have suggested that the merger 
fraction may peak at some critical redshift. However, un- 
til now there have been very few observational studies for 
z > 1 to constrain the redshift of this possible peak, and 
hence the detailed physics of the galaxy merging at high 
redshift is still uncertain. 

The standard hierarchical formation model has been 
modified to incorporate the feedback effects of a central 
SMBH (eg. Silk & Rees 1998). Hydrodynamical simula- 
tions have suggested that as galaxies merge, they evolve 
into a starburst system, and a substantial amount of gas 
is eventually funneled onto their SMBHs. This triggers 
a visible AGN typically 1-2 Gyr after the onset of the 
merger (Springel, Di Matteo, & Hernquist 2005). While 
previous authors have discussed an evolutionary link be- 
tween mergers and starburst galaxies (eg. Larson & Tins- 
ley 1978; Sanders & Mirabel 1996; Bell et al. 2005), this 
"duty cycle" of merging, followed by a starburst, followed 
by an AGN, has only recently been suggested (Hopkins et 
al. 2006). This formation scheme implies that we should 
observe similarities in the redshift evolution of merging 
galaxies, their star formation history (eg. Madau et al. 
1996), the mass assembly in galaxies (eg. Dickinson et al. 
2003), their AGN activity (eg. Ueda et al. 2003; Hasinger, 
Miyaji, & Schmidt 2005; Schneider et al. 2005), and re- 
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suiting galaxy morphologies (eg. Driver et al. 1995; Abra- 
ham et al. 1996). Since these other phenomena have been 
well studied at z>l, we desire corresponding measure- 
ments of galaxy merging at comparable redshifts, which 
is the focus of this work. 

While there may be some convergence on the the- 
oretical front, there is little consensus regarding the 
observational connections between merging and active 
galaxies. In general, AGN host galaxies do not ex- 
hibit strong merger-type morphologies (eg. Bahcall et al. 
1997; Cohen et al. 2006; Guyon, Sanders, & Stockton 
2006) or show enhanced image asymmetry (Grogin et al. 
2005). Conversely, Canalizo et al. (2007) find evidence 
for recent merger activity in the quasar host galaxy, 
MC2 1635-1-119, which was previously classified as an 
undisturbed elliptical galaxy. Some of this disagreement 
may be resolved in the context of the above "duty cy- 
cle," wherein the morphological indications of the recent 
merger have subsided by the time the AGN is observa- 
tionally identified. This debate will surely continue ow- 
ing, in part, to the difficulties interpreting morphological 
studies (eg. de Propris et al. 2007). 

This work is organized as follows; In § 2 we describe the 
dataset and the mass estimates, § 3 we outline our merger 
selection, § 4 we measure the merger fraction and num- 
ber density, § 5 we compare this study to published work, 
and conclude with a general discussion of the result in the 
context of galaxy evolution in § 6. We adopt the three- 
year results from the Wilkinson Microwave Anisotropy 
Probe (Spergel et al. 2007), where ^^0 = 0.24, i7A = 0.76, 
and Ho = 100h km s~^ Mpc~^ with /i = 0.73. All magni- 
tudes quoted herein are in the AB system (Oke & Gunn 
1983). 

2. DATA 

The deep fields observed by the Hubble Space Telescope 
(HST) provide an unprecedented view of the distant uni- 
verse (z > 1), and are therefore the ideal datasets for 
studying the high redshift evolution of the galaxy merger 
rate. The GRism ACS Program for Extragalactic Sci- 
ence (GRAPES; Pirzkal et al. 2004) has provided low- 
resolution {R ~ 100) spectroscopy of 1400 objects with 
z' < 27.2 mag in the Hubble Ultra Deep Field (HUDF; 
Beckwith et al. 2006). From these data, Ryan et al. 
(2007) computed spectrophotometric redshifts for 1308 
galaxies at z < 5. When the Balmer or 4000 A break 
occurs in these optical spectra (for 0.5 < z < 1.5), the 
standard deviation in (zphot — •Zspcc)/(1 + Zgpcc) is 0.037, 
and it is only 0.050 for z > 1.5, which make these among 
the most accurate photometric redshifts computed using 
the standard maximum-likelihood technique (Bolzonella, 
Miralles, & Pello 2000) to date. 

Since the primary aim of this work is to measure the 
high redshift evolution of the galaxy merger rate, we re- 
quire a set of criteria to select our merger candidates. 
While this has traditionally been performed using the 
rest-frame luminosity, we will select our systems based 
on their mass. Ideally, we desire their total gravitational 
mass, since it is expected to drive galaxy merging (eg. 
Khochfar & Burkert 2001). However, since this mass is 
exceedingly difficult to determine for a large sample of 
faint galaxies, we will adopt the stellar mass as our se- 
lection quantity. 

To determine the required stellar masses of our sample, 
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Fig. 1. — The merger candidate selection. We show the com- 
puted stellar mass as a function of the spectrophotometric redshift. 
The Ughtly sliadcd region indicates the volumc-liniitcd, mass selec- 
tion region for tlic primary galaxies, while the darker region shows 
the possible selection of companion galaxies. In total, ~ 95% of 
the galaxies with z' < 27.2 mag identified in the HUDF by Beck- 
with ct al. (2006) are represented here (sec Figure 2) with 230 of 
these meeting the selection criteria discussed in § .3. The various 
lines show the mass of a maximally old, passively-evolving stellar 
population with an e-folding time in the star formation rate given 
in the lower-right corner. 

we fit the combination of GRAPES spectra and broad- 
band observations in the B-, J-, and i?-bands with a 

grid of stellar population synthesis models (Bruzual & 
Chariot 2003, BC03 hereafter). The masses are deter- 
mined by minimizing the reduced between the spec- 
trophotometry and these models, following the method 
of Papovich, Dickinson, & Ferguson (2001). By parame- 
terizing the star formation history of these galaxies as 
5'(t) ~ ^'oe~*/'^, our model grid consists of three pa- 
rameters: the extinction in the restframe V-band, the 
e-folding time in the star formation history, and stellar 
population age. If we assume a Salpeter initial mass 
function (Salpeter 1955) and solar metallicity, then the 
parameter space consists of ~4 x 10^ independent mod- 
els. By varying the values of the grid points, we estimate 
the uncertainty on the stellar masses to be < 15%. In Fig- 
ure 1, we show the stellar mass as a function of redshift 
for these galaxies. The different tracks show the mass of 
a maximally old stellar population for the various values 
of T listed in the lower-right corner. The shaded region 
shows the mass selection for the primary galaxies, which 
is ~95% complete for z' <27.2 mag. The darker area in- 
dicates the selection region for the lower mass companion 
galaxies to be discussed below. 

3. IDENTIFICATION OF MERGING GALAXIES 

Since we have reasonable redshift information for the 

entire sample of galaxies, it is natural that we adopt the 
dynamically close pair method (eg. Carlberg, Pritchet, 
& Infante 1994). However, since our redshift uncertainty 
is roughly a factor of 10 larger than that of a tradi- 
tional spectroscopic survey, our approach will be some- 
what akin to the more common pair count analysis (eg. 
Zepf & Koo 1989). This technique has several key virtues 
not present in other approaches: it has simple and well 
understood statistics to compute (Abraham 1999) and 
does not rely on calibration datasets (de Propris et al. 
2007). Patton et al. (2000, P2000 hereafter) demonstrate 
that the galaxy merger rate ("T^mg) can be computed from 
the galaxy merger or pair fraction (/) and the galaxy 



High Redshift Galaxy Mergers 



3 



; HUDF Umll I i 

: GRAPES Umit i • : 

3 0.6 H I H 

f OA - I I -- 

0.2 H , H 

O.O E- ;.. V ^ , : V • • 

20 22 24 26 28 30 32 34 
i' (mag) 

Fig. 2. — The i'-band number counts. Since the GRAPES objects 
were selected from the HUDF to have z' < 27.2 mag, it is not 
surprising that > 95% complete at these flux levels. In the lower 
panel, we show the ratio of the number counts of the GRAPES 
survey to all objects in the HUDF. 

number density (n) as TZmg = fnr^g , where Tmg is the 
typical merger timescale. 

3.1. Pair Criteria 

Using the redshift catalog of Ryan et al. (2007), we 
require the following four criteria to be satisfied to iden- 
tify a galaxy pair as a merger system: (1) at least one 
neighboring galaxy <20 kpc away to ensure that 
close galaxy pairs will merge (eg. de Propris et al. 2005; 
Bell et al. 2006) on a timescale <0.5 Gyr (eg. Khochfar 
& Burkert 2001; Patton et al. 2002); (2) we require the 
mass ratio of primary galaxy to the putative companion 
to be < 4 : 1 , which is consistent with a flux ratio between 
the galaxies of <1.5 mag; (3) we require the stellar mass 
of the primary galaxy to be Af* > 10^" M0; and (4) wc 
require the rcdshifts of the primary and the companion 
to agree within the 95% confidence interval of the spec- 
trophotometric redshifts of Ryan et al. (2007). 

3.2. Accounting for the Flux and Mass Limit 

While the GRAPES survey is complete to very faint 
flux levels (sec Figure 2), we may expect it to be volume- 
limited at high redshift and that no corrections to the 
pair statistics are needed. To investigate this notion, we 
show the absolute i?-band magnitude as a function of 
the spectrophotometric redshift in Figure 3 where the 
solid line shows our flux limit of z' ~ 27.2 mag for a 
flat spectrum in following Driver (1999). While this 
figure shows the expected volume-limited behavior (eg. 
de Propris et al. 2005), there is a more subtle question 
to address since this is a mass-selected sample: Is there 
a population of galaxies which are fainter than our flux 
limit and more massive than our mass limit? If such a 
population exists, then our mass selection would be in- 
complete in these galaxies (Dickinson et al. 2003). At the 
redshifts of interest, such objects must have a relatively 
high stellar mass-to-light ratio {'^*,b ^2 Mq Lq^) and 
therefore be very red in the HST images. 

To find these dark, red galaxies at 1.5<z<2.4, Daddi 
et al. (2004) have developed a technique based on the ob- 
served {B — z') and {z' — K) colors (the BzK method). 
When this is applied to the HUDF, Daddi et al. (2005) 
identify only seven old, passively-evolving galaxies with 
K <22) mag. They argue that these systems must have 
been formed in a quick burst (r < 0.3 Gyr), since they 
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Fig. 3. — The B-band absolute magnitude shown for compar- 
ison and not used in the merger selection. The solid line shows 
the contour of constant apparent magnitude of z' ~ 27.2 mag 
following Driver (1999). The majority of our galaxies have ages 
<500 Myr, which implies their median stellar mass-to-light ratio 
is T,,B !^0.4 Mq Lq^. Therefore, the mass limit of M, > lO^o Mq 
roughly corresponds to Mg < —20.5 mag, which is comparable to 
the luminosity limit imposed by various authors (eg. Patton et al. 
1997; Le Fevre et al. 2000) at lower redshifts. 

emit little rcstframe, ultraviolet fiux. While these objects 
meet the above mass-to-light criteria, they are brighter 
than our flux limit and are included in our redshift sur- 
vey. Therefore, to estimate the number of galaxies satis- 
fying this mass-to-light criterion and are fainter than our 
flux limit, we must perform a similar analysis on deeper 
data. Since there is no sufficiently deep K-h&nd imag- 
ing in the HUDF, we must use the NICMOS imaging 
in the i7-band (Thompson et al. 2005) and assume that 
these galaxies typically have {H — i^T) ~0.8 mag (Forster 
Schreiber et al. 2004). In Figure 4, we show the corre- 
sponding color-color space for the BzH selection for the 
objects in the NICMOS imaging with fluxes from Coe et 
al. (2006). The shaded polygon shows the BzH selection 
region defined by {z' -H)-{B-z') < -0.2-(H-K) mag 
and {z' — H) > 2.5 — {H — K) mag. The open squares rep- 
resent all objects satisfying these BzH criteria while the 
open circles show the five objects of Daddi et al. (2005) 
which are in this portion of the HUDF. As in Figure 1, 
the various tracks show BC03 stellar population mod- 
els as indicated in the lower- left. Since we recover 4/5 
objects from Daddi et al. (2005), we conclude that this 
is a viable method to find the maximally old, passively- 
evolving galaxies which satisfy the mass-to-light require- 
ment. There are only 9 galaxies (or 1% of all objects 
with z'<28.5 mag) which are expected to be burst-type 
systems (or a t<0.3 Gyr; Daddi et al. 2005). However, 
the abundance of galaxies with masses and colors com- 
parable to the maximally old, passively-evolving models 
with r>0.3 Gyr in Figure 1 and Figure 4 suggests that 
corrections to the pair statistics are necessary. 

To correct the observed pair counts, we will adopt the 
method of P2000, and for brevity only highlight our mod- 
ifications for the mass-selection and give the essential de- 
tails. In the P2000 method, one computes the maximal 
A"-corrcx;tion from a set of galaxy SED templates to derive 
the low luminosity limit when selecting objects based on 
absolute magnitudes (see Figures 1 and 2 from Patton et 
al. 2002). However, since we select candidates based on 
their stellar mass, we must determine the corresponding 
maximal mass limit. While Patton et al. (2002) adopt 
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Fig. 4. — The BzH selection for passively evolving, burst-type 
galaxies. By assuming {H — K) ~ 0.8 mag (Forstcr Sclirciber et al. 
2004), we perform a similar analysis of Daddi et al. (2005) on the 
much deeper NICMOS imaging in the HUDF to estimate the range 
of star formation histories applicable to these redshifts. The shaded 
region represents the BzH selection discussed in § 3.2 and, as in 
Figure 1, the different lines represent maximally old, passively- 
evolving stellar populations, denoted by the e-folding time in their 
star formation rate in the lower-right. The "buUseye" marks the 
birth of the populations at Zform^oo- As discussed in § 3.2, maxi- 
mally old galaxies with t < 0.3 Gyr are exceedingly rare and would 
appear in the shaded, BzH polygon. Therefore, when correcting 
our pair counts following the P2000 method, we adopt BC03 mod- 
els appropriate for our sample {r>0.3 Gyr). 

the empirical templates of Coleman, Wu, & Weedman 
(1980), we choose the stellar population synthesis models 
of BC03. As in § 2, we parameterize the star formation 
rate as ^'(t;T) oc e"*/"^ and restrict t > 0.3 Gyr, since 
objects with r<0.3 Gyr are extremely rare (as discussed 
above). Since the aim is to generate a population of max- 
imal stellar mass, we set t = T{z), where T{z) is the age 
of the Universe at rcdshift z. 

We define the stellar mass limit to be the more massive 
between the fiducial limit (M,,o) and the maximal mass 
model (M*(z;t)): 



Mi™(z) = max [M,,o, M,(z; t)]. 



(1) 



With this mass selection, we compute the selection func- 
tion to correct the counts of the flux-limited survey as 



Sn{z) = 



Mi'"(z) 
M..0 



'^{M^,z)dM^ 



J^^^^<^>{M,,z)dM. 



(2) 



where ^{M^,z) is the galaxy mass function (not to be 
confused with the luminosity function). Eisner, Fcul- 
ner, & Hopp (2007) show that the galaxy mass function 
can be well modeled as the Schcchter function (Schechter 
1976) with M^^IO"-"^ Mq and a ~ -1.36 for z<5. Fi- 
nally, the optimal weights to apply to the pair statis- 
tics will be wpf(z) oc Sn{z)^^, where the various re- 
maining terms (such as boundary effects and spectro- 
scopic completeness) arc computed following Patton et 
al. (2002). Using these optimal weights, the number of 
close companions for the ith primary galaxy is computed 
as Nc^i = WN^lzj), where j runs over the compan- 
ions which satisfy the above merger criteria. Finally, the 
mean number of close companions per primary galaxy is 
found as 

^^^ Z^nM)N (3) 
following the notation of Patton et al. (2002). 



3.3. Error Budget 

While it may be appropriate to assume Poisson uncer- 
tainties on the pair counts, we require a more detailed 
analysis to incorporate the effects of the redshift uncer- 
tainties. Despite having redshifts for each component 
of the merging pair, the relatively large uncertainties on 
the proper distance {6dp < 500 Mpc) can lead to false 
positives. To estimate the magnitude of this effect, we 
randomize the redshifts and masses in such a manner to 
preserve the selection discussed above and recompute the 
pair counts. In 10^ realizations of the redshift catalog, we 
estimate that <1.5 pairs per redshift bin may arise from 
superpositions identified by erroneous redshifts. This ap- 
proach allows us to bootstrap the uncertainties on the 
pair counts to include the effects of Poisson shot noise, 
chance superpositions, and redshift uncertainties. 

Since the HUDF covers ~ 3.'3 x 3.'3, the uncertainties 
resulting from cosmic variance must also be addressed. 
By assuming the correlation function of the galaxies is 
^(r) oc (r/ro)^'*', then the relative cosmic variance in the 
number counts is given as = 1.45(V/rQ)~'''/^ for cell 
volume V, correlation length ro, and 7=1.8 (Somerville 
et al. 2004). For comparable sampled volumes, Beckwith 
et al. (2006) estimate the cosmic variance to be cr^ ~ 0.17 
for the HUDF cell geometry. 

4. REDSHIFT EVOLUTION 

In Figure 5, we show the resulting galaxy merger frac- 
tion as a function of {1 + z) as filled circles with un- 
certainties that reflect the contributions from chance su- 
perpositions, cosmic variance, redshift uncertainty, and 
Poisson shot noise. Since these galaxies are typically very 
young (an age < 500 Myr as inferred from the SED fit- 
ting described in § 2), their median stellar mass-to-light 
ratio in the restframe S-band is T*^b ~ 0.4 Mq Lg^, 
and therefore the mass limit of M*,o > 10^° Mq roughly 
corresponds to Mb ^ —20.5 mag. Since this luminosity 
limit is similar to that imposed by Patton et al. (1997) 
and Le Fevre et al. (2000), we show those measure- 
ments for comparison as open triangles and diamonds, 
respectively. We show a standard parameterization of 
/(2) =0.02 X (1 + z)'^ as a dashed hue. While this model 
can fit the low redshift {z < 1) observations, it fails to 
reproduce our higher redshift {z>l) data. Therefore, we 
also show the empirical formed proposed by Conselice 
(2006) of 

f{z) = ail + z)'e'^^'+^\ (4) 

as a solid line. Using a standard least-squares technique, 
we find a = 0.44 ± 0.04, b = 7.0± 0.6, and c = 3.1 ± 0.4 
for xt = 0.28. This model has a maximum at Zf,ac = 
(5/c) — 1 = 1.3 ± 0.4, when accounting for the correlated 
uncertainties on b and c. In the inset, we show the 68%, 
95%, and 99% confidence intervals in the b-c plane for 
the fit to these data. The peak in the fraction is likely 
a consequence of the strict mass limit and mass ratio 
imposed, and not a result of flux incompleteness (see 
Figure 1). Our observations can rule out the standard 
power-law model which has a reduced x^ = 7.8. 

While the merger fraction may provide excellent infor- 
mation on galaxy evolution, the galaxy merger number 
density is more relevant when comparing to other extra- 
galactic observations. From the merger fraction (/), we 
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Fig. 5. — The observed galaxy merger fraction. The results from 
this work are indicated by filled circles, while the observations of 
Patton et al. (1997) and Le Fevre et al. (2000) arc shown by open 
triangles and diamonds, respectively. The dashed line shows the 
common form of f{z) = 0.02 X (1+^)^ , which fails to fit our high red- 
shift observations. However, the model f(z) = a{l + z)^e~'^^^'^'\ 
shown as the solid line, can reproduce all the observations rea- 
sonably well (xf, = 0.28). This peak in the merger fraction at 
^frac — 1.3it0.4 may arise from the strict mass limit (M, >10"' Mq) 
and mass ratio (< 4 ; 1) imposed on our sample and is not due to 
flux incompleteness given the unparalleled depth of this survey (see 
Figure 1) and the corrections applied in § 3.2. 
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Galaxy Merger Results 



2 










(10-4 Mpc-3) 


0.75±0.25 


0.09±0.03 


6.31±2.34 


1.25d=0.25 


0.11±0.04 


3.56±1.52 


2.00±0.50 


0.09±0.06 


2.04±1.39 



t Computed from the merger fraction by equation (5). 

compute the merger number density as 



Jz—Az 

where N{z) is the completeness corrected redshift distri- 
bution from the Ryan et al. (2007) catalog. Then the 
merger rate is TZmg = J^mgTmg (P2000). In Figure 6, we 
show the merger number density for the same data as 
discussed in Figure 5. As a solid line, we show empir- 
ical model of equation (4), converted to a number den- 
sity using the redshift distribution of equation (23) in 
Bcnftez (2000). The merger number density peaks at a 
Zdens — 1-2, which is a slightly lower redshift than the 
peak in the merger fraction. Finally, we tabulate our 
measured galaxy merger fraction and number density in 
Table 1. 

5. COMPARISON TO PREVIOUS WORK 

For low to moderate redshifts (z < 1), there have 
been numerous measurements of the galaxy pair fraction. 
While many studies disagree on the precise values of the 
pair fraction, they generally agree that the fraction rises 
from ~ 2% at 2 ~ to ~ 15% at 2 ~ 1 as a power-law in 
(1-1-2;). Where our observations overlap with previous 
work, we find a merger fraction of 15±5% at z^l, which 
is in good agreement with published values (eg. Le Fevre 
et al. 2000; Lin et al. 2004; Kartaltepe et al. 2007). 



Fig. 6. — The observed galaxy merger number density. As in 
Figure 5, the results from this work are indicated by filled circles, 
and the observations of Patton et al. (1997) and Le Fevre ct al. 
(2000) are shown as open triangles and diamonds, respectively. The 
dotted line shows the total density of AGN with Lx > 10**^ erg s" ^ 
selected in the soft (0.5—2 keV) X-ray band (Hasinger, Miyaji, & 
Schmidt 2005). The number density of merging galaxies peak at 
•2dons — 1-2 while the AGN number density peaks at zagn — 0.65. 
This difference in redshift corresponds to a cosmic time difference 
of ~ 2.5 Gyr. A comparable time delay between the merging and 
the AGN phase has been suggested by simulations (Springcl, Di 
Matteo, & Hcrnquist 2005). Moreover, the ratio of the densities 
at their peaks would suggest that as much as ~ 20% of the major 
merging of massive galaxies will evolve into X-ray AGN with Lx > 
10*2 erg s"^ (eg. Daddi et al. 2007). 

At higher redshift {z>l), there are no other studies of 
merging using the pair fraction, due to the difficulty in 
acquiring sufficiently deep redshift datascts. However, if 
the redshift of the primary galaxy can be obtained, then 
a morphological approach can be a viable method for de- 
termining the merger status. Conselice et al. (2003) has 
shown a roughly power-law form to z~3 using the CAS 
(concentration, asymmetry, dumpiness) morphological 
system (Bershady, Jangren, & Conselice 2000; Conselice 
2003). When these measurements are coupled with rest- 
frame colors, the refined merger fraction is substantially 
higher than our observations (Conselice 2006). However, 
de Propris et al. (2007) discuss some of the hmitations of 
this system and caution that merger studies using image 
asymmetry require careful interpretation. 

Based on the matched Two Micron All Sky Sur- 
vey (2MASS)/Two-Degrcc Field Galaxy Redshift Sur- 
vey (2dFGRS; Cole et al. 2001), Xu, Sun, & He (2004) 
identify 38 merging systems with similar close pair cri- 
teria at 0.01 < 2 < 0.07. By adopting a standard stellar 
mass-to-light ratio in the restframe /Cg-band i^*,Ks = 
1.32 M0 Lq""^), these galaxies arc selected to have > 
-|^q10.25 While they restricted their work to a smaller 
range of mass ratios for their systems {5Ks < 1 mag), they 
find a pair fraction of 1.70 ±0.32%, consistent with other 
values at these redshifts (eg. Patton et al. 1997). 

6. DISCUSSION 

Using the deepest imaging and spectroscopic data 

available, we show that the fraction and number den- 
sity of major mergers between massive galaxies peak at 
relatively low redshifts. These peaks are likely not an ar- 
tifact of the data nor the measurement process, given the 
careful attention to the mass limits, spectrophotometric 
redshift accuracy, accounting of measurement uncertain- 
ties, unprecedented flux limit, and the good agreement 
with published data in the redshift range of overlap. 
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Since the high redshift portion of our observed merger 
rate dramatically differs from what is found by numerous 
authors, it is interesting to explore some of the conse- 
quences which arise from this difference. In particular, 
one implication of our merger rate concerns the fraction 
of galaxies which have undergone a major merger, or 
the remnant fraction /rem- By assuming a merger frac- 
tion of f{z) ~ 0.011 X (1 -I- z)^'^ and a merger timescale 
r„,g = 0.5 Gyr, Patton et al. (2002) find that /rem -15% 
of galaxies with — 21 < Mb < — 18 mag have experienced 
a major merger since This is in contrast with the 

semi-analytic models of De Lucia et al. (2006), who argue 
that ~ 50% of stars which comprise elliptical galaxies as- 
semble into a single dark halo between < z < 1 by major 
merging. When we apply this analysis^ to our slightly 
larger and steeper merger rate, we find that f^^m — 42% 
of galaxies with > 10"'^" M0 have been involved in a 
major merger since z — 1. Moreover, since our merger 
rate peaks at some critical redshift, the fraction of mas- 
sive galaxies that have undergone a major merger tends 
to a limiting value at 2; —» 00 of /ji™ ~ 62%, whereas a 
strict power-law form increases to 100%. This suggests 
that (1 — /'em) — 38% of massive galaxies may have never 
undergone a major merger in their lifetime, but does not 
preclude a series of minor mergers. 

The magnitude of physical processes, which are driven 
by galaxy merging, should be proportional to the merg- 
ing galaxy number density. Many phenomena have been 
linked to galaxy merging, most notably, star formation, 
AGN activity, and galaxy morphologies. The cosmic 
star formation rate density, as traced by radio, ultra- 
violet, far infrared, and/or various emission line lumi- 
nosities has been measured by numerous authors at z < 7 
(Lilly et al. 1996; Madau et al. 1996; Connolly et al. 
1997; Steidel et al. 1999; Haarsma et al. 2000; Brinch- 
mann et al. 2004; Wyder et al. 2005; Bouwens et al. 
2006; Hopkins & Beacom 2006). While there is rea- 
sonable consensus that the star formation rate density 
increases from z ~ to z ~ 1, the high redshift portion 
can differ greatly among the many authors. Nonethe- 
less, the peak in our merger number density is consistent 
with the peak in the cosmic star formation rate density 
at z~ 1—1.5. Moreover, the peak star formation rate den- 
sity is ~10"i M© yr~i Mpc"-'' (eg. Hopkins et al. 2006), 
which would require the typical merger system to have a 
star formation rate of ~100 Mq yr~^. This rate is con- 
siderably higher than what is expected for major mergers 
(Bell et al. 2005), which suggests that major mergers arc 
not the only source of star formation at these redshifts. 
However, minor merging is expected to be more frequent, 
and may be responsible for triggering the remaining por- 

^ Our remnant fraction is found from equation (18) of Patton et 
al. (2002). 



tion of the cosmic star formation density. 

Many authors have addressed an evolutionary link be- 
tween merging galaxies and AGN (Carlberg 1990; Silk 
& Rees 1998; Kauffmann & Haehnelt 2000; Croton et 
al. 2006; Hopkins et al. 2006), and generally agree that 
major merging can drive AGN activity. Therefore, it is 
natural to compare our merger number density to the 
observed AGN density. The density of type-1 AGN se- 
lected in soft (0.5-2.0 keV) X-ray band with luminosi- 
ties of Lx = 10^^ — 10'*^ erg s~^ from Hasinger, Miyaji, 
& Schmidt (2005) is shown in Figure 6 as a dotted line. 
This sample provides the roughly instantaneous AGN ac- 
tivity unobscured through any dust, and is similar to 
the results in the hard X-rays (Ucda et al. 2003). This 
AGN density peaks at 2;agn — 0.65, or a lookback time of 
~5.9 Gyr, whereas our observed merger number density 
peaks at J^icns — 1-2, or ~ 8.5 Gyr in lookback time. If 
a major merger of two galaxies triggers an AGN, then 
we expect on average ~2.6 Gyr between the time when 
the galaxies are observed as a pair on the sky, and when 
the remnant would be observed as an AGN in the soft 
X-rays. Springel, Di Matteo, & Hernquist (2005) have 
suggested that major merging can trigger the fueling of 
the central black hole after the merger onset, and will be 
observed as a visible AGN ~ 1—2 Gyr later. According to 
these simulations, we expect the merger number density 
to peak ^2.5 Gyr earlier in cosmic time than the AGN 
density which is similar to our observations. 

These measurements have exploited the excellent red- 
shift accuracy, spatial resolution, and depth of ACS grism 
observations. The Wide Field Camera 3 for HST will 
provide combined ultraviolet prism, infrared grism. and 
numerous ultraviolet, optical, and near-infrared broad- 
band observations. With such a rich dataset, redshifts 
of comparable accuracy from 2 < 2; < 6 should become 
standard. This will be ideal for refining this work, and 
determining the role of various effects, such as galaxy 
environments and/or cosmic variance at z:<6. 

We are very grateful for the excellent comments and 

suggestions of the anonymous Referee which have im- 
proved this paper. We additionally thank Rolf Janscn, 
Norman Grogin, Nimish Hathi. Sangeeta Malhotra, 
James Rhoads, and Nor Pirzkal for a number of helpful 
discussions. RER acknowledges support from the ASU 
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